Coupled phantom field in loop quantum cosmology 
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A model of phantom scalar field dark energy under exponential potential coupling to barotropic 
dark matter fluid in loop quantum cosmology is addressed here. We derive a closed-autonomous 
system for cosmological dynamics in this scenario. The expansion in loop quantum universe has 
a bounce even in presence of the phantom field. The greater decaying from dark matter to dark 
phantom energy results in greater energy storing in the phantom field. This results in further turning 
point of the field. Greater coupling also delays bouncing time. In the case of phantom decaying, 
oscillation in phantom density makes small oscillation in the increasing matter density. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

There has recently been evidence of present acceler- 
ating expansion of the universe from cosmic microwave 
background (CMB) anisotropics, large scale galaxy sur- 
veys and type la supernovae [3, [2, Q . Dark energy (DE) 
in form of either cosmological constant or scalar field 
matter is a candidate answer to the acceleration expan- 
sion which could not be explained in the regime of stan- 
dard big bang cosmology DE possesses equation of 
state p = wp with w < —1/3 enabling it to give repulsive 
gravity and therefore accelerate the universe. Combina- 
tion of observational data analysis of CMB, Hubble Space 
Telescope, type la Supernovae and 2dF datasets allows 
constant w value between -1.38 and -0.82 at the 95 % 
of confident level [f|. Meanwhile, assuming flat universe, 
the analysis result, —1.06 < w < —0.90 has been re- 
ported by [6j using WMAP three-year results combined 
with Supernova Legacy Survey (SNLS) data. Without 
assumption of flat universe, mean value of w is -1.06 
(within a range of -1.14 to -0.93). Most recent data (flat 
geometry assumption) from ESSENCE Supernova Sur- 
vey la combined with SuperNova Legacy Survey la gives 
a constraint of w — — 1.07±0.09 7]. Observations above 
show a possibility that a fluid with w < — 1 could be al- 
lowed in the universe Q . This type of cosmological fluid 
is dubbed phantom. Conventionally Phantom behavior 
arises from having negative kinetic energy term. 

Dynamical properties of the phantom field in the stan- 
dard FRW cosmology were studied before. However the 
scenario encounters singularity problems at late time 
[l(J. While investigation of phantom in standard 
cosmological model is still ongoing, there is an alterna- 
tive approach in order to resolve the singularity prob- 
lem by considering phantom field evolving in Loop Quan- 
tum Cosmology (LQC) background instead of standard 
general relativistic background [ll|, • L°°P Quantum 



Gravity-LQG is a non-perturbative type of quantization 
of gravity and is background-independent [HI, [lij] • LQG 
provides cosmological background evolution for LQC. An 
effect from loop quantum modification gives an extra cor- 
rection term —p 2 /p\ c into the standard Friedmann equa- 
tion [l5|,[l6|,[l7j. Problem for standard cosmology in dom- 
ination of phantom field is that it leads to singularity, so 
called the Big Rip [18]. The —p 2 /p\ c term, when dom- 
inant at late time, causes bouncing of ex pan sion hence 
solving Big Rip singularity problem [l4l Il9l. [2fj| . Re- 
cently, a general dynamics of scalar field including phan- 
tom scalar field coupled to barotropic fluid has been in- 
vestigated in standard cosmological background. In this 
scenario, the scaling solution of the coupled phantom 
field is always unstable and it can not yield the observed 
value ~ 0.7 [2l|. Indeed there should be other effects 
from loop quantum correction to the Friedmann equa- 
tion. Moreover when including potential term in scalar 
field density, the quantum modification must be included 
[22l ]. Although, the Friedmann background is valid only 
in absence of field potential, however, investigation of 
a phantom field evolving under a potential could reveal 
some interesting features of the model. In this letter, 
we investigate a case of coupled phantom field in LQC 
background in alternative to the standard relativistic cos- 
mology case. In Section UU we introduce framework of 
cosmological equations before considering dynamical au- 
tonomous equations in Section IIIIl We show some nu- 
merical results in Section ITVl where the coupling strength 
is adjusted and compared. Conclusion and comments are 
in Section W\ 
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II. COSMOLOGICAL EQUATIONS 
A. Loop quantum cosmology 

The effective ffat universe Friedmann equation from 
LQC is given as fTEL [l7l] . 



H 2 = 
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1 - 



P_ 
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(1) 



where H is Hubble constant, Mp is reduced Planck mass, 
p is density of cosmic fluid, pic = V3/(16ttC 3 G 2 ?i). The 
parameter £ is Barbero-Immirzi dimensionless parameter 
and G is the Newton's gravitational constant. 



B. Phantom scalar field 

Nature of the phantom field can be extracted from ac- 
tion, 
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Energy density and pressure are given by 
1 



p^ = --r + V(&, 



and 



with equation of state, 



W<f, 
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When the field is slowly rolling, the approximate value 
of w is -1. As long as the approximation, <j) 2 ~ or the 
condition, <fi 2 < 2V holds, w is always less than -1. In 
our scenario, the universe contains two fluid components. 
These are barotropic fluid with equation of state p m — 
p-mWra and phantom scalar field fluid. The total density 
is then p = p m + p,j, which governs total dynamics of the 
universe. 



C. Coupled phantom scalar field 

Here we consider both components coupling to each 
other. Fluid equations for coupled scalar fields proposed 
by [iH assuming flat standard FRW universe are 



f3F(l 
3if(H 
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(6) 
(7) 



These fluid equations contain a constant coupling be- 
tween dark matter (the barotropic fluid) and dark energy 



(the phantom scalar field) as in [2J|. Eqs. ([6]) and ([7]) 
can also be assumed as conservation equations of fluids 
in the LQC. Total action for matter and phantom scalar 
field is lH 



S 



d 4 a 



M 2 



+ Sm(0). (8) 



Assuming scaling solution of the dark energy in standard 
cosmology, the pressure Lagrangian density is written as 



v{X,<t>) = -X- cexp(-A<^/M|), 



(9) 



where X is the kinetic term, —g ab d a (j)db(f>/2 of the La- 
grangian density © and ©. The second term on 
the right of Eq. Q is exponential potential, V{4>) — 
cexp(— \(f)/Mp) which gives scaling solution for canoni- 
cal and phantom ordinary scalar field in standard general 
relativistic cosmology when steepness of the potential, A 
is fine tuned as 



A = Q 



1 + W m - Q^jlVm - W,p) 



(10) 



The steepness (fTU]l is, in standard cosmological circum- 
stance, constant in the scaling regime due to constancy 
of and £1$ [HI, HI]- However, in LQC case, there has 
been a report recently that the scaling solution does not 
exist for phantom field evolving in LQC Therefore 
our spirit to consider constant A is the same as in [25[ 
not a motivation from scaling solution as in [23|. The 
exponential potential is also originated from fundamen- 
tal physics theories such as higher-order gravity [26( or 
higher dimensional gravity [27| |. 



III. COSMOLOGICAL DYNAMICS 

Time derivative of the effective LQC Friedmann equa- 
tion LQC CQ) is 



2M| 
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In above equations we define new variable 

S = 4>. (14) 

The coupled fluid equations © and (f?]) are re-expressed 
in term of S as 

dV 

S = -3HS+ — + Qp m , (15) 
<±<p 

p m = ~3H(l + w m )p m + Qp m S . (16) 



FIG. 1: Phase portrait of S(t) versus (j>(t) for Q = 
—0.4, —0.1, 0.1 and 0.4 from left to right. All trajectories have 
the same initial conditions 5(0) = 0.5 and 0(0) = 0. 
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FIG. 2: Scale factor plotted versus time for Q = 
-0.4, —0.1, 0.1 and 0.4 (from bottom to top). 



The Eqs. {T3]), (0, |T5]) and ([16]) form a closed au- 
tonomous set of four equations. The variables here are 
p m , S, 4> and H. The autonomous set recovers standard 
general relativistic cosmology in the limit p\ c — ► oo. The 
general relativistic limit affects only the equation involv- 
ing H. From the above autonomous set, one can do a 
qualitative analysis with numerical integration similar to 
phase plane presented in different situation [28| . Another 
approach of analysis is to consider a quantitative analysis 

m 



FIG. 3: Phantom field density plotted versus time for Q = 
—0.1. The other values of Q also yield bouncing and oscilla- 
tion. 

IV. NUMERICAL SOLUTIONS 

Here we present some numerical solution for a positive 
and negative coupling between the phantom field and 
barotropic fluid. The solutions presented here are phys- 
ically valid solutions corresponding to Class II solutions 
characterized in [ill ]. For non- minimally coupled scalar 
field in Einstein frame [3fj| , the coupling Q lies in a range 
-1/V6 < Q < 1/V6. Here we set Q = -0.4,-0.1,0.1 
and 0.4 which lie in the range. Effect of the coupling can 
be seen from Eqs. © and ([7]). Negative Q enhances de- 
cay rate of scalar field to matter while giving higher mat- 
ter creation rate. On the other hand, positive Q yields 
opposite result. Greater magnitude of Q < gives higher 
decay rate of the field to matter. Greater magnitude of 
Q > will result in higher production rate of phantom 
field from matter. 



A. Phase portrait 

The greater Q value results in greater value of the field 
turning point (see ^-intercept in Fig. [1]). The kinetic 
term S(t) turns negative at the turning points corre- 
sponding to the field rolling down and then halting before 
rolling up the hill of exponential potential. When Q is 
greater, the field can fall down further, hence gaining 
more total energy. The result agrees with the prediction 
of Eqs. © and 0. 

B. Scale factor 

From Fig. [21 the bounce in scale factor occurs later for 
greater Q value of which the phantom field production 
rate is higher. The field has more phantom energy to 
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accelerate the universe in counteracting the effect of loop 
quantum (the bounce). For less positive Q, the phantom 
production rate is smaller, and for negative Q, the phan- 
tom decays. Therefore it has less energy for accelerating 
the expansion in counteracting with the loop quantum 
effect. This makes the bounce occurs sooner. 



C. Energy density 

Time evolutions of energy density of the matter and 
the phantom field are presented in Figs. [3]and[31 If Q > 
0, the matter decays to phantom. This reduces density 
of matter. While for Q < 0, the matter gains its density 
from decaying of phantom field. In Fig. [3] there is a 
bounce of phantom density before undergoing oscillation. 
For a non-coupled case, it has recently been reported that 
the phantom density also undergoes oscillation [Tl|. As 
seen in Figs, m and 03 the oscillation in phantom density 
of the phantom decay case (Q < 0) affects in oscillation in 
matter density while for the case of matter decay (Q > 0), 
the matter density is reduced for stronger coupling. The 
oscillation in the phantom density comes from oscillation 
of the kinetic term </>, i.e. S(t) as shown in [ill ]. 

V. CONCLUSION AND COMMENTS 

In this letter, we have derived an autonomous system 
of a loop quantum cosmological equations in presence of 
phantom scalar field coupling to barotropic matter fluid. 
We choose constant coupling Q between matter and the 
phantom field to positive and negative values and check 
numerically the effect of Q values on (1) phase portrait, 
(2) scale factor and (3) energy density of phantom field 
and matter. We found that field value tends to roll up 
the hill of potential due to phantom nature. With greater 
Q, the field can fall down on the potential further. This 
increases total energy of the field. For canonical scalar 
field either standard or phantom, LQC yields a bounce. 
The bounce is useful since it is able to avoid Big Bang 
singularity in the early universe. Here our numerical re- 
sult shows a bouncing in scale factor at late time. This is 
a Type I singularity avoidance even in presence of phan- 
tom energy. The greater coupling results in more and 
more phantom density. Greater phantom effect therefore 
delays the bounce, which is LQC effect, to later time. In 
the case of matter decay to phantom (Q > 0), oscilla- 
tion in phantom energy density does not affect matter 
density. On the other hand, when Q < 0, phantom de- 
cays to matter, oscillation in phantom density results in 
oscillation in the increasing matter density. 

This work considers only the effects of sign and mag- 
nitude of the coupling constant to qualitative dynamics 
and evolution of the system. Studies of field dependent 
effects of coupling Q(<fi) in some scalar-tensor theory of 
gravity and investigation of an evolution of effective equa- 
tion of state could also yield further interesting features 
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FIG. 4: Matter density plotted versus time for Q = 
—0.4, —0.1, 0.1 and 0.4 (from top to bottom). 
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FIG. 5: Zoom-in portion of Fig. [4] The phantom field decays 
to matter at highest rate for Q = —0.4 (top line). Oscilla- 
tion in matter density due to oscillation in the phantom field 
density is seen clearly here. 



of the model. Quantitative dynamical analysis of the 
model under different types of potential is also motivated 
for future work. Frequency function of the oscillation in 
scale factor and phantom density are still unknown in 
coupled case. It looks like that the oscillation frequency 
tends to increase. This could lead to infinite frequency 
of oscillation which is another new singularity. 
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